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In the paper " Superconducting decay length in a ferromagnetic metal" by Gusakova, Kupriyanov 
and Golubov [Pis'ma v ZhETF 83, 487 (2006); cond-mat/0605137 , the authors claim that they 
solved the linearized Eilenberger equation in the ferromagnetic region of an S/F heterostructure at 
arbitrary mean free path. In this comment we show that the solution suggested by the authors is 
not correct and explain details of the exact solution found by us in an earlier work several years ago 
(Ref.[2]). 



I. INTRODUCTION 



In a recent paper 0] Gusakova, Kupriyanov and Golubov studied theoretically the decay length over which the 
condensate function / decreases in the ferromagnetic region of a S/F structure (S and F denote a superconductor 
and ferromagnet). In order to determine the decaying length they have solve the linearized Eilenberger equation 
representing the solution in the form 

f(x) = C(9)exp(-x/0, (1) 

where 9 is an angle between the direction of the momentum and the a;-axis (the x-axis is perpendicular to the S/F 
interface). In this solution the function C(0) depends only on 8 and the decay length £ is assumed to be independent 
of this angle. Substituting this Ansatz into the linearized Eilenberger equation the authors obtain expressions for 
C{9) and the decay length £. They analyze the dependence of the decay length on the exchange field h and the mean 
free path for an arbitrary magnitude of the product hr (we set the Planck constant equal to 1). 
In this Comment we would like to point out that: 

1. The solution of the form JQ) for the linearized Eilenberger equation is not correct and may serve only for rough 
estimates for the decay length £ in a particular limit. 

2. The linearized Eilenberger equation can be solved exactly for an arbitrary value of product hr. The exact solution 
has been obtained in Ref. for arbitrary values of hr (contrary to the statement of the authors of Ref. that the 
solution was obtained in [2| in the clean limit only) and used later in subsequent theoretical [3|, |^| and experimental 
studies [ip. 

In Ref. [21 the solution of the linearized Eilenberger equation for the Josephson S/F/S junction is presented explicitly. 
In order to avoid a literal repetition of that result we use the method of Ref. Q , to obtain the exact solution for a 
simpler S/F structure considered in thus directly comparing the result with "the solution" of Eq. Q from that 
work. From a general formula we obtain simple expressions for f(x) in the limiting cases of small and large values of 
the product hr (the diffusive and quasi-ballistic cases correspondingly). To make the analysis more transparent, we 
neglect the spin-orbit interaction. 

The solution of the linearized Eilenberger equation f(x) can be represented as a sum of a symmetric s(x) and 
antisymmetric a(x) parts: f(x) — s(x) + a(x). The functions f{x), s(x) and a(x) are matrices in the particle-hole 
space. The antisymmetric part a(x) is related to the symmetric one via the formula 

a = -sgnw (fil/Ku,) t 3 8 x s (2) 
and the symmetric part s(x) in the ferromagnetic region (x > 0) obeys the equation 0] 
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H I d x J - kJ = -k u (s) (3) 

where / = vpT is the mean free path, /i = cos#, n u = (1 + 2|o; m |r) — sgno;2zft,r, u) m = nT(2m + 1) is the Matsubara 
frequency. The angle brackets mean the angle averaging 

(s) =(1/2) £ dtf(j*) (4) 
Eqs J2J-ISI are complemented by the boundary condition Q 

a | x =o+= 7(^)sgnw (f 3 f s J , (5) 

where = T(/x)/4 and T(/z) is the transmission coefficient. The matrix f s — 1T2A/ y/uj^ + A 2 = i^fs is the 

matrix condensate function of the superconductor (the BCS-function). The boundary condition (J5J is valid in the 
case of a small transmission coefficient T(/i) when the linearization of the Eilcnberger equation is possible. 

We can formally continue the solution of Eq. ((3J) into the region of negative x in a symmetric way taking into 
account the boundary condition (J5J. Therefore Eq. Q may be written in the form 

» 2 l 2 d* x s - k Is = -k u [(s) + 2 1 (^f s S(x)} (6) 

where the function s(x) is defined as s(x) — s{x)iT2. The boundary condition (JjjJ is satisfied automatically due to the 
last term on the right-hand side. For the Fourier transform Sf. = J dxs(x) exp(—ikx) we get from Eq.@) 

«*M = J^tt^)) + 2 7(^)^/ s } , (7) 

where M{n) = (klfi) 2 + k 2 ,. 

Averaging Sfe(/i) over \x and substituting the result into Eq.J7J), we obtain finally 

2k w £/ s k u 7(/i)/i 

Sfc(Al) = w { rMi7M) ( ir 1 + 7( ^ } (8) 

The inverse Fourier transformation 

/dk 
—s k (fj,)exp(ikx) (9) 

determines the exact solution for Eq.© at arbitrary mean free path. It is obvious that the function s(x,/^) cannot 
be represented as a product of two functions, such that one of them would depend only on /i and the other one only 
on x. 

The average (1/M) can easily be calculated 

1 kl 

(1/M) = ^T tan "- (10) 

Let us consider limiting cases, where one can obtain analytical expressions for the condensate function s(x,/z). 
a) diffusive limit: hr << 1. 

In this case, characteristic wave vectors k are much smaller than X/l. Expanding the function on the right-hand 
side in Ea. ljTU|l in powers of kl, we obtain 1 — n u (l/M) sw l 2 (k 2 + k%)/3. The behavior of s(x,fj,) at large x (larger 
than the mean free path I) is determined by the pole in the first term in Eq.©. Calculating the residue at this pole, 
we get 



3 



s(x,fx) « -r -^{y(fx)iJ,)exp(-k h x) (11) 
k h l 

where fc^ = 2[|u m | — sgnui m ih\/ D : D = lv/3 is the diffusion coefficient. This result can be easily obtained from the 
Usadel equation. One can see that s(x, //) is small and the linearization of the Eilenberger equation is possible provided 
the condition 3v / 3(7(/i)A 1 ) /V~hT << 1 is fulfilled. Only in this case the characteristic length of the condensate decay 
(at large a;) does not depend on the angle 6. At x of the order of the mean free path the solution of Eq.© depends 
on 8. Thus, a solution of the form Eq. is only valid in the diffusive limit and for large x. However, in the case of 
strong ferromagnets or not too short mean free path the solution depends on angles at any x, as we will show, 
b) quasi-ballistic case: hr >> 1. 

In this case the parameter | k w |w hr is large and therefore the first term in the curly brackets in Eq. (JHJ) is small. 
Calculating the integral in Eq.©, we get 



s(x,[i) « 7(/i)/ s exp(-/c w a;/^), ^>0 (12) 

We see that in this case the decay length £p which is equal to £p = min{u/x/T, fil} is angle dependent. The 
condensate function oscillates in space with a small (compared to £f) period £ osc = l/i/ (2irhr). It is clear that the 
condensate function can not be presented in the form 

In a general case of arbitrary product hr (hr w 1) the spatial dependence of the Cooper pair wave function / can 
be found by numerical integration of Eas. H7l8(l . 
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